Mixing structures in the Mediterranean Sea from Finite-Size Lyapunov Exponents 
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We characterize horizontal mixing and transport structures in the surface circulation of the 
Mediterranean Sea, as obtained from a primitive equation circulation model. We calculate the 
Finite Size Lyapunov Exponents (FSLEs) of the velocity data set, which gives a direct measure of 
the local stirring. By proper election of the FSLE parameters, we focus on the mesoscale structures, 
locating a number of vortices embedded in an intricate network of high-stretching lines. These lines 
control transport in the system. At the edge of the vortices, a dense tangle of line intersections ap- 
pears, identifying strong mixing. The spatial distribution of FSLEs, averaged over one year, allows 
to classify areas in the Mediterranean basin according by their mixing activity. The space average of 
FSLEs on selected geographical regions gives a measure for quantifying and comparing the mixing 
seasonal variability. 
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Introduction 

Horizontal transport and mixing processes are central 
to the study of the physical, chemical, and biological dy- 
namics of the ocean. Correct understanding and precise 
modelling of them are relevant from a theoretical view- 
point and crucial for a range of practical issues, such 
as plankton dynamics or the fate of pollutant spills. In 
this regard, the last few years hav e seen the appearance 
of interesting new developments [Mariano et al.. 2002j 
on the Lagrangian description of transport and mix- 
ing phenomena, many of them coming from the area 
of nonlinear dynamics. Such approaches do not aim 
at predicting individual tracer trajectories, but at lo- 
cating spatial structures that are known from dynam- 
ical systems theory to act as tem plates for the whole 
flow [Ottino. 19 89. Wi aains. 1992j| . Such structures (at- 
tractors, saddles, manifolds, ...) have been used since 
many years ago for classifying the evolution of trajec- 
tories in abstract dynamical systems. However, when 
put in a fluid dynamics context, they gain a new and 
direct physical meaning, corresponding for instance to 
avenues and barriers to transport, vortex boundaries, or 
lines of strong stretching. In particular, chaotic motions 
such as the ones occurring in the turbulent ocean are 
characterized by complex intersection of stretching and 
contracting manifolds at the so called hyperbolic points: 
regions of fluid initially compact in the proximity of these 
points become elongated along the stretching directions 
and then folded, leading to the typical filamental and con- 
voluted structures that are common in satellite pictures 
of water temperature or chlorophyll, as well as in labora- 
tory experiments. From the point of view of transport, 
such regions are characterized by a strong mixing: tra- 
jectories of initially close particles are quickly separated 
along the stretching directions, and fluid of different ori- 
gins is inserted in between. 

Until recently, the power of these novel Lagrangian ap- 
proaches has been mainly relegated to mathematical sys- 
tems or simplified workbench models, since the required 



detailed knowledge of the velocity field was not read- 
ily available in real geophysical situations. However, in 
the last decades the situation has dramatically changed, 
with a rapidly increasi ng; amount of data available 
from Lagrangi an drifters \Mariano et al. . 2002j . satellite 



measurements 



Ealvcrn. 20Q PJ, and especially from de- 
tailed compute r models \Haidvoael and Beckmann. 19981 
\Dietrich. 1997| . This has paved the way to a grow- 
ing number of geophysical applications, among which 
we men tion studies of the Loo p Current in the gulf of 
Mexico Kuznetsov et al.. 2002j . or the characterization 
of dispersion properties in the Adriatic Sea from drifter 
experiments y|. The methodology is not restricted to 
ocean circulation, but it is also being developed in atmo- 
spheric dynamics Q, and mantle convection 0]. Indeed, 
new Lagrangian approaches are very appealing for geo- 
physical applications, since they can be used as a tool for 
automatically extracting transport structures underlying 
raw Eulerian velocity data. 

In this Letter we characterize mixing strength at the 
mesoscale in different areas of the Mediterranean Sea by 
means of a Lagrangian technique, the Finite Size Lya- 
punov Exponents method. The technique also identifies 
dynamical objects that organize the transport, and rele- 
vant coherent structures. 



Numerical Data 

We produce velocity data from the DieCAST ocean 
model (Dietrich for Center Air Sea Technology), adapted 
to the Mediterranean basin with a horizontal resolu- 
tion of 1/8 deg rees (approximately 10 km) and 30 verti- 
cal levels (see Fernandez et al.. 2004j for details). The 
DieCAST model is a z-level primitive equation model 
based on the hydrostatic, incompressible, and rigid lid 
approximations. The model has been integrated for 
20 years being forced by yearly repeating monthly cli- 
matological atmospheric forcing. Using such climato- 
logical forcing the model reproduces well the general 
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surface circulation and many of the important features 
of the observed annual cycle of the Mediterranean Sea 
Fernandez et al.. 2004| . Due to the adequate horizon- 



tal resolution and the numerical characteristics of the 
model, basically high ord er numerics an d low numerical 
and physical dissipation Dietrich. 1997| . the numerical 
simulations reproduce a great number of mesoscale struc- 
tures, that are required for the present study. 

In this work we analyze the daily model output of ve- 
locity field obtained from the last three years of simu- 
lation and corresponding to the second vertical horizon- 
tal layer (in this way we avoid the strong dependence on 
wind stress of the first model layer). This layer has a ver- 
tical extent of 11.60 to, centered at a depth of 16 to. The 
two-dimensional horizontal velocity field on this layer is 
not exactly incompressible, but it is very close to this 
situation since typical vertical velocities in the ocean are 
four orders of magnitude smaller than horizontal ones 
(10 -5 vs 10 _1 to/s). Thus, points at which fluid particles 
diverge horizontally along particular directions receive 
fluid along other directions, so that stretching is linked to 
mixing. On the other hand, within the Finite-Size Lya- 
punov exponent (FSLE) method used in this Paper and 
presented in the next sections to estimate transport at 
the mesoscale, fluid particle trajectories need to be inte- 
grated only while they undergo horizontal displacements 
of the order of 100 km, i.e. during 1-10 days (see Section 
of Results). Estimating an effective or average vertical 
velocity for this time as the spatial average of the ver- 
tical velocity field in horizontal regions of that size, one 
obtains an effective vertical velocity of 0.1 — 0.7 m/day. 
Thus, during the time of integration, most particles do 
not leave the horizontal layer considered. In consequence, 
restricting the study to horizontal motion on a single 
model layer is a good description of the full transport pro- 
cesses for the space and time scales relevant to mesoscale 
processes within the FSLE approach. 



Finite Size Lyapunov Exponents (FSLEs) 

A common way to quantify the stretching by advec- 
tion is by means of the standard Lyapunov exponents. 
They arc defined as the exponential rate of separation, 
averaged over infinite time, of fluid parcels initially sep- 
arated infinitesimally. In realistic situations (such as the 
case of the Mediterranean Sea where boundaries at finite 
distance strongly influence the circulation) the infinite- 
time limit in the definition makes the Lyapunov expo- 
nent a quantity of limited practical use. Recently, the 
Finite Size Lyapunov Exponent (FSLE) has been intro- 



duced Aurell et al. l99li\Artale et al. 1997j in order to 



ing and visualizing Lagrangian structures (e.g., transport 
barriers or vortex boundaries) Q. Here we will focus 
mainly in the second use, but we will also introduce mea- 
sures of dispersion and mixing based on the Lagrangian 
structures detected. 

The FSLE technique appears to be ideally suited 
for oceanographic applications, being the mathematical 
analogous of a floater experiment: a set of tracers, with 
some initial mutual distances, are followed in time as 
they are transported by integrating the velocity field (we 
use a bilinear interpolation to assign velocities to points 
that are not model grid points; thus the velocity field is 
effectively smooth at scales below 1/8 of degree). The 
FSLE is inversely proportional to the time at which two 
tracers reach a prescribed separation. More precisely, 
A(x, t,do,5f), the FSLE at position x and time t, is com- 
puted from the time r it takes for a trajectory starting 
at time t at a distance So from x to reach a separation 
Sf from the reference trajectory that started at x: 



1 Sf 

X(x,t,S ,Sf) = - log-/-. 

t o 



(1) 



study non-asymptotic dispersion processes, which is par- 
ticularly appropriate to analyze transport in closed areas. 
FSLEs have been used for two complementary goals: for 
characterizing dispersion processes U and for detect- 



In order to characterize the strongest separation (and 
the fastest convergence along the complementary direc- 
tion) , A is selected as the maximum among the four values 
obtained when the initial separation So is chosen along 
four orthogonal directions. 

The FSLE depends on the choice of two length scales: 
the initial separation So and the final one Sf. Here, we 
are interested in the spatial distribution of FSLEs, and 
thus we calculate them at points x located on a grid of 
spacing Ax. In this case, a simple argument shows the 
convenience of using a value of So close to the intergrid 
spacing Ax: If one chooses Sq much smaller than Ax, all 
the points of a stretching manifold laying further than Sq 
from any grid point are not tested, and thus the method 
gives only a rather discontinuous sampling of the struc- 
ture. On the other hand, if So is much larger than Ax, the 
same stretching manifold is detected ( "smeared" ) on sev- 
eral grid points, with a loss in spatial resolution. Since we 
are interested in mesoscale structures, the other length, 
Sf, will be chosen as Sf = 1 degree, i.e., separations of 
about 110 Km. In this way the FSLE represents the 
inverse time scale for mixing up fluid parcels between 
length scales S and Sf. 

Results 

The spatial distribution of FSLEs for a particular day 
of the simulation (10th of June of the first year of the 
data set) can be seen in Fig.^ Typical values are in the 
order of 0.1 — 0.6 days -1 , corresponding to mixing times 
for mesoscale distances of 1.7 — 10 days. As observed 
in previous works, maximum values of the distribution 
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organize in lines [2], |3j that provide good approximations 
to repelling material lines (which are in turn stable man- 
ifolds of hyperbolic moving points) . These lines orga- 
nize the transport processes in the basin. Spatial struc- 
tures ranging from the small scales to the ones typical 
of mesoscale vortices are clearly identified. Computing 
such picture for every day of the year and taking then 
the time average, one can obtain a map of regions in the 
Mediterranean with different mixing activity (fig- ■ As 
expected, the Southern part of the basin appears more 
active, especially close to the North African coast. In 
fig. |3] we show more in detail the FSLEs of the area in 
the small box of fig. In the core of the eddies one has 
low values of the FSLEs (i.e., low dispersion rates); on the 
contrary, the largest values of the FSLEs can be found in 
the outer part of the eddies, where the stretching of the 
fluid parcels is particularly important. 

Note that in some regions of the vortex cores, chaotic 
tangles are still observed as local maxima of the FSLE 
distribution. These maxima are, however, not strong. 
In fact, even if the stretching is locally very high, the 
requirement for two points to diverge for more than 5f = 
110km gives a low value of A to such finer structures, 
since such distance is bigger than the size of the vortex 
that acts as a containing barrier most of the time. This is 
an example of a useful property of the FSLE technique: 
it allows to restrict the analysis to the structures relevant 
for transport among selected lengthscales only. 

Additional sets of coherent structures and organizing 
lines can be obtained by computing the FSLEs from tra- 
jectory integration backwards in time. Maxima in the 
new distribution identify lines of maximum compression, 
approximating attracting material lines or unstable man- 
ifolds of hyperbolic moving points . Since stable and 
unstable manifolds cannot be crossed by particle trajec- 
tories, such lines strongly constrain and determine fluid 
motion. 

Calculating in this way the FSLEs in a region of strong 
mixing, we unveil the tangle of stretching and compress- 
ing lines in which vortices are embedded (fig. |3J . These 
lines also define the directions of transport. Lobes arising 
from intersections of stretching and compressing lines at 
a vortex edge indicate where transport in and from the 
vortex takes place, whereas tangencies among them pro- 
vide barriers to transport. In Fig.0]some intersections of 
stretching and compressing lines are indicated as black 
dots. These identify Lagrangian hyperbolic points (and 
their motion define hyperbolic trajectories). Such points 
correspond to areas with strong mixing activity: fluid 
is advected here along a compression line and then dis- 
persed away along the stretching line. 

This dynamical picture suggests a quantitative mea- 
sure of mixing in a prescribed area A: one can define de 
quantity M±(t) = (-^/A+A_) A where A+ and A_ are the 
FSLEs in the forward and in the backwards time direc- 
tion, and the average is the spatial average over the area 



A. This quantity is large only where hyperbolic points 
are present. The time dependence of this quantity when 
the area A is the whole Mediterranean is shown in fig. [5] 
characterizing the seasonal variations of mixing. Max- 
imum values, of the order of 0.13 days -1 , are attained 
in winter. Because of the approximate incompressible 
character of the horizontal flow, the temporal variations 
of forward and backward FSLEs are strongly correlated, 
and one expects that the same information can be ob- 
tained from just one of the FSLEs. Thus one can define a 
simpler measure of mixing in an area as M + (t) = {\+} A - 
We show in fig.|5]that, as expected, it contains essentially 
the same information as M±, but at variance with it, it 
could in principle be measured from floater experiments. 
It is thus a more convenient characterization of mixing 
strength. As a further example, we compare (fig. EJ) 
the temporal behavior of M + (t) in two regions (the areas 
North and South of the Balearic islands, corresponding 
to the boxes in fig. |5J) where we expect (from fig. [2J to 
see a very different mixing activity. The higher activ- 
ity in the Southern part, where the Algerian current is 
present, is confirmed. In addition, seasonal fluctuations 
are smaller in the Northern part. 



Conclusions 

The FSLEs provide a direct method for computing si- 
multaneously the mixing activity and the coherent struc- 
tures that control transport at a given scale. Interest- 
ingly, the method is based on the evolution of the relative 
separation between two passive tracers, and thus numer- 
ical results directly suggest floater-based experiments for 
verification. In this work we have analyzed with this 
method horizontal velocity data from a Mediterranean 
computer simulation. Different mixing behavior between 
geographical regions, and at different seasons, is readily 
characterized. Finally, we point out that the strong hor- 
izontal stirring associated with hyperbolic points should 
also have important biolo gical conseq uences. In the di- 
rection of recent works \Mo,riin. 2 003] . it would be inter- 
esting to compare the Lagrangian structures presented 
here with productivity, patchiness, or other measures ob- 
tained from biological distributions. 
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FIG. 1: FSLE spatial distribution for the whole Mediter- 
ranean on a specific day (10th of June of the first simulation 
year). Ax = So = 0.02 degrees. Mesoscale structures and 
vortices can be clearly seen. Units for FSLEs are day -1 . A 
zoom of the indicated box is presented in fig. |H] 
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FIG. 2: Time average (for the first simulation year) of the 
FSLEs in the whole Mediterranean basin. Geographical re- 
gions of different mixing activity appear. Colors as in fig. 
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FIG. 3: Enlarged plot of the small box shown in fig.Q 
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FIG. 4: FSLEs calculated from forward (displayed as positive 
values) and backwards (displayed as negative values) integra- 
tions in time. A region with strong mixing appears organized 
by a tangle of stretching and compressing manifolds. Such 
lines organize the flow. The black dots indicate some of the 
hyperbolic points that are located at the intersections of the 
lines. 
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FIG. 5: Temporal evolution of the mixing measures M+(t) 
(circles), and of M±(t) (squares) for the whole Mediterranean 
Sea during one simulation year. They display a very similar 
behavior, with maximum values in winter. 
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FIG. 6: M+(t) during one year for the Algerian current 
(squares), which corresponds to the area delimited by the 
southern box in fig. |5] and the north of the Balearic islands 
(circles), which is the northern box in fig. [5] 



